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ABSTRACT: We show how to simulate U(l) gauge fields coupled to three-dimensional 
quantum gravity and then examine the phase diagram of this system. Quenched mean field 
theory suggests that a transition separates confined and deconfined phases (for the gauge 
matter) in both the negative curvature phase and the positive curvature phase of the 
quantum gravity, but numerical simulations find no evidence for such transitions. 



Introduction 



In two dimensions, dynamically triangulated random surfaces have proven to be a 
useful discretization of quantum gravity. An analogous approach in three dimensions has 
been formulated [1,2]. However, this system can only be studied in the grand canonical 
ensemble because no set of ergodic moves have been found for fixed volume. The action 
can be written as 

S = aN - /3N 3 (1) 

where Nq is the number of nodes and N% is the the number of tetrahedra (i.e. the volume). 
(3 is therefore the cosmological constant while a acts like Newton's constant, since at 
fixed volume varying the number of nodes varies the curvature. In two dimensions, the 
topological classification of manifolds allows one to prove that at fixed Euler character 
there is an exponential bound on the number of triangulations as a function of the number 
of triangles. This implies that a chemical potential can control the volume. In three 
dimensions there is no corresponding classification and therefore no corresponding proof of 
an exponential bound. The naive result is factorial growth which would mean no partition 
function could be defined. It was a significant discovery that a chemical potential (/?) does 
control the volume in three dimensions [3,4]. The volume diverges at a particular critical 
value of (3. Given that the thermodynamic limit can be taken, the next issue is the nature 
of the phase diagram. In particular, is there a second order phase transition where a non- 
trivial continuum limit could be taken? It turns out that there are two phases separated 
by a transition at a*. For smaller values of a the system is in a negative curvature phase 
where the coordination number of the nodes diverges. For larger values of alpha, the 
system is in a positive curvature phase where the number of sites is proportional to the 
volume. Unfortunately, the transition between these two phases, where one might have 
hoped for a continuum limit with zero curvature, appears to be first order so that there is 
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no interesting continuum limit. 

This doesn't necessarily spell doom, since an expanded phase diagram (with more 
coupling constants) might still contain an appropriate transition. A simple and interesting 
way to expand the phase diagram is to couple matter to the system. Recently, studies of 
Ising matter coupled to three-dimensional quantum gravity have been completed [5-7] . If 
the matter is placed on the nodes, the phase diagram consists of the gravitational transition 
and an Ising transition. The two cross at zero Ising coupling so that the Ising matter does 
not change the nature of the gravitational transition. If a dual formulation is employed in 
which the matter fields are placed on the elementary volumes, the gravitational transition 
is also unchanged. 

In this paper, we couple gauge fields to three-dimensional quantum gravity. (Our 
simplicial lattices correspond to triangulations of S^). Aside from the possibility of finding 
a non-trivial continuum limit, this system is interesting because it couples gauge fields 
and quantum gravity. In two-dimensions, the analogous system is topological — quantum 
gravity has no interesting effect [8]. Three dimensions is thus the lowest dimension in 
which coupling these two theories gives a non-trivial result. In what follows, we restrict 
ourselves to the U(l) gauge theory. 



A quenched mean field theory proved useful in the case of Ising matter coupled to 
three-dimensional quantum gravity and it proves useful for gauge matter as well. The 
action for U(l) gauge fields coupled to three-dimensional quantum gravity is 



A 

where the sum is over all triangles and is the sum of gauge angles about the given 
triangle. A simple approximation to the phase diagram can be obtained by doing mean 



Mean Field Theory 




(2) 
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field theory for the gauge variables while fixing the volume and the number of nodes. This 
allows us to consider the influence of the background metric on the matter, but neglects 
the influence of the matter on the metric. The simplest mean field theory [9] for gauge 
theories begins with a trial action 

S H = Hj2 cos6 ( 3 ) 

links 

for which the free energy per link is simple 

u(H) = lnI (H) (4) 

where Iq is a modified Bessel function. Adding and subtracting the trial action to the 
gauge part of the full action gives 

Z = Tre S = Z(H) < exp(A ^ cos <p A - H ^ cos 9) > H 

A links 

> Z(H)exp < A^cos^a - H ^ cos/9 > H (5) 
A links 

where Z is the partition function for the gauge part of the full action, Z(H) is the partition 
function for the trial action and the subscript H denotes that the trial action is used as 
the weight in the expectation values. The resulting free energy per unit volume is 

F(H, A) = (uNi - Hu'Ni + X(ufN 2 )/N 3 (6) 

(the prime denotes the derivative with respect to H, N\ is the number of links, and 
N 2 is the number of triangles). This must be maximized with respect to the field H. 
Maximizing gives a first order transition at A* = 0.794 in the negative curvature phase 
and A* = (1 + 1/c) * 0.794 in the positive curvature phase, where c = N%/Nq as the 
volume is taken large. (A more sophisticated gauge invariant approach leads to the same 
results). This result suggests the thick line in the phase diagram in figure 1. Note that 
mean field theory cannot be taken too seriously and is given here just to guide the later 
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numerical calculations. In three dimensions on a fixed lattice, for instance, the theory can 
be rewritten as a Coulomb gas that is in an ionized phase for all temperatures. Mean field 
theory mistakes a crossover for a phase transition. 

Placement of the gravitational transition in the phase diagram is guided by contact 
with the pure gravity theory that can be made in the limits of large and small A. Note 
that (using N\ = Nq + JV3) 

-2n 



links 

so that 



J] Qf "tWj = (2tt)^ = e^ ln27r = e (^o+iV3)ln2 7 r 



(7) 



Z(a,p,X = 0) = Z purcgrav (a + ln27r,/3-ln27r) (8) 

implying that the gravitational transition is shifted to a* ~ 1.91 — 2.16 [4]. For very large 
A the theory also approaches a pure gravity theory. In this limit, (since N2 = 2./V3) 

A cos(0 A ) ~ XN 2 = 2AiV 3 (9) 

A 

implying 

Z(a, f3, A) ~ -^pure grav(«, P — 2A) (A such that cos <j)/± ~ 1) (10) 

This final equation suggests that the gravitational transition is only weakly dependent on 
A as indicated by the dotted line in figure 1. 



Detailed Balance 



The task now is to numerically check figure 1 and to determine the order of the phase 
transitions suggested by mean field theory. In order to simulate this theory, we must derive 
the detailed balance relations in the grand-canonical ensemble. In three dimensions, any 
triangulation can be reached from any other triangulation through a series of local moves 
chosen from a set of four possibilities. These moves are labelled by the number of tetrahedra 
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before and after the move. The (1,4) move (figure 2) corresponds to the insertion of a node 
into a tetrahedron and the (4,1) move is its inverse. The (2,3) move (figure 3) replaces a 
triangle separating two tetrahedra with a new link at the boundary of three tetrahedra. 
The (3,2) move is its inverse. In contrast to the situation with Ising matter, where only 
the (1,4) detailed balance relations were modified by the addition of matter, here all the 
relations are modified. First, consider the (2,3) move. Here, it is necessary to assign a 
value for the gauge field on a potential new link. Label the state with two tetrahedra A 
and the state with the new link 6>, where 9 is an angle between and 2tt. Then consider 
attempting a 2 — > 3 move with 50% probability and a 3 — > 2 move with 50% probability. In 
the 2^3 case choose the triangle by picking a tetrahedron at random and then choosing 
one of its four sides. Note that there are two ways of picking the same triangle with this 
procedure. In the 3^2 case, the link is chosen from a list of all possible links with 
coordination number 3. If 1% is the number of such links, the resulting detailed balance 
equation is 

The acceptance rate can be increased by integrating over the possible values of 9 so that 
it is first decided whether to update the triangulation or not and then later the gauge 
field values are determined. Integrating over all possible values of 9 and assuming that 
P{9 — > A) can be chosen independent of 9 gives 

-^rP(A -> B) = P(B -> A) e S0 - SA d9 (12) 
2^3 Jo 

where B denotes the state with a new link regardless of its value. Define 

hz = / e Se ~ SA dff (13) 
JO 

(which will turn out to be a gauge invariant function of the loops on the surface of the 
cluster of the two or three tetrahedra under consideration) and denote the ratio 

#23 = h/2N s (14) 
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then detailed balance is satisfied by 



P[A - B) = (16) 
P{B - A) = (16) 

Now consider detailed balance for the (1,4) move. In this case, attempt a 1 — > 4 move 
or a 4 — > 1 move each with a 50% probability. In the former case, just pick a tetrahedron 
at random. In the latter case, pick from a list of sites known to have coordination number 
four. Detailed balance requires 

lw/ Ap{A ~> e) = lr/ ep{e ~> A) (17) 

where A represents the state with no node, O represents the four gauge degrees of freedom 
living on the four new links, and Z4 represents the number of nodes with coordination 
number four. Again, define 

I u = J e S@ ~ SA (18) 
where now integration is over four angles, and 

R U = h/N 3 (19) 

then detailed balance is satisfied by 

P(A - B > - R^ki (20 > 
p{B ^ A)= whi (21) 

where now B represents the presence of a node without specifying the four new gauge 
degrees of freedom associated with it. These update probabilities tell us whether to put 
in or take out a link and whether to put in or take out a node. They do not indicate what 
the gauge degrees of freedom should be. This is taken care of by a subsequent Metropolis 
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or heat bath step. One caveat should be added. The ratio in -R14 is interpreted differently 
depending on whether a 1 — > 4 move is being contemplated or a 4 — > 1 move is being 
contemplated. In the former case, ./V3 is the number of tetrahedra to begin with while I4 is 
the number of sites with coordination number 4 if the move is accepted. In the latter case, 
I4 is the number of sites with coordination number 4 to begin with while N3 is the number 
of tetrahedra there will be if the move is accepted (i.e. the current number of tetrahedra 
minus 3). Similar considerations are required in the case of the (2,3) move. 

At first sight, these equations suggest that a simulation of gauge fields coupled to 
quantum gravity in this grand-canonical ensemble would be impractical. An integration 
needs to be done every time there is a local update. We will see that the integral J23 is 
not really very troublesome. The integral /14 is much more demanding computationally, 
but it turns out to be within the realm of practicality. 

Consider first the action difference for the former integral (the (2,3) move). 



where 1, 2, and 3 label the new triangles associated with the new link and labels the 
old triangle that has been eliminated. The quantities t{ are the cosines of the sum of the 
gauge angles about the specified loop. The needed gauge integral is therefore 



JO 

If the links, ctj, involved in the integral are labeled as in figure 3 and X{ is defined such 
that 



(3 + A(*i + 1 2 + 1 3 - t ) 



(22) 




(23) 



x\ = a\ + ct2, X2 = «3 + «4, £3 = «5 + olq 



(24) 



then 



^1+^2 + ^3 = 003(6* + x\) + cos(6 ) + X2) + cos(6> + £3) 



(25) 



r cos(6 ) + uj) 



(26) 
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where u> is independent of 9 and 



3 \ 2 / 3 - 



r 2 



c ° s x i\ + [Yl sin x * ) ( 27 ) 

u=l / \i=l 



Finally, 

7 23 = 2tt 7 (Ar) e^^ (28) 

Is this result gauge invariant? The values of the gauge fields enter only through r. It is 
easy to see that r is gauge invariant. Change variables to 

z\=x\- X2, Z2 = x\ - X3, ,23 = X2 - £3 (29) 

then 

r = 3 + 2(cos z\ + cos Z2 + cos z%) (30) 
Since each of the Z{ is a gauge loop, r is manifestly gauge invariant. 

il4 is more complicated. Here there are four gauge variables to be integrated over, 
labeled as in figure 2. The action difference is 

~~ = ct — 3(3 + X cos Pi (31) 

i=\ 

where the pi are the new plaquettes. In terms of the gauge variables 

pi = 6>i + 2 + xi (32) 

P2 = 01 + % " ^5 (33) 

P3 = 01 - 04 " ^3 (34) 

P4 = 02 ~ 03 - £4 (35) 

P5 = 02 + 04 " ^2 (36) 

P6 = 03 + 04 " ^6 (37) 
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The action ratio can now be rewritten 

z=6 / oo \ 
e S e -S A = e «-3/J "Q J- I k .(X)e ikiPi \ (38) 

i=l \ki=— oo / 
The exponential can be rewritten as i times 

+ k 2 + h) + e 2 (h + k 4 + k 5 ) + e 3 (k 2 - fc 4 + fc 6 ) + ^(-^ + h + h) 

+x\ki — x<^k 2 — x^ks — X4/C4 — x 2 k§ — xqUq (39) 
Integration over 6\, 6 2: and #3 gives a factor 

(2tt) 3 5(h + k 2 + fc 3 ) 5(A;i + A; 4 + fc 5 ) 5(^2 - A; 4 + fc 6 ) (40) 

Notice that these three 5-functions automatically make the coefficient of #4 zero. This is 
a consequence of gauge invariance under transformations at the center node. Integration 
over #4 just gives another 2n. The 5-f unctions give 

h = ~h - k 2 , &5 = -k\ - £4, kQ = /C4 - /C2 (41) 

so that the xi terms become 

k\{x\ + x 2 + X3) + k 2 (-x§ + X2, + xq) + k4(-X4 + x 2 - xq) (42) 

Each coefficient of the k^ is a gauge invariant loop on the surface of the tetrahedron. If 

the coefficients of the k{ are relabeled as z\, z 2 , and z% respectively then the final result is 

00 00 00 

iu=(2^-^ y: e e w a ) 

mi—— 00 m 2 =— 00 m 3 =— 00 

/ m2 (A)/ m3 (A)/ mi+m2 (A)/ mi+m3 (A)/ m3 _ m2 (A)e^^^+-^2+m 3 z 3 ) (43) 

Invariance under sign changes in mi, m 2 , and 7713 make it clear that the right hand side 
is real — the exponent can be replaced with a cosine. (Computations are faster if the 
expression is rewritten to take advantage of this fact) . 
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Results 



One parameter that has to be tuned is the ratio, r, of update sweeps of the gauge 
field to the update sweeps of the metric. In two dimensions, when bosons are coupled to 
quantum gravity, the expectation value of the bosonic action can be calculated exactly due 
to scale invariance. This can be used to check whether the bosons are in equilibrium and 
then the equilibrium of the gravity sector can be checked by varying the above ratio and 
checking to see that the results are stable. Here, there is no analytical check of either the 
matter or the gravity sector. We can only check various expectation values as a function 
of the ratio and look for stability. For a fixed lattice (that has not been warmed up) one 
update pass on the gauge fields is insufficient. We suspect that on a dynamical lattice, 
where the connectivity is constantly changing, that more than one gauge update is required 
per update sweep through the metric. Consider the expectation values of the plaquette, 
the specific heat, and the number of nodes, as well as the value of (3 required to maintain 
a given volume. As r is varied from 1 to 16, the values of all of these quantities change, 
but reach a plateau near r = 8. This is illustrated in fig. 4 with a plot of the expectation 
value of the number of nodes versus r. The expectation value of the plaquette has a similar 
behavior, varying from 0.62 to 0.67 as r varies from 1 to 16, but varying by less than a 
part per thousand as r is increased from 8 to 16. r = 8 is the value we used. A number of 
runs at r = 1 were also done. The results differed numerically from the r = 8 results, but 
they implied the same conclusions described below. 

The phase diagram suggested by the quenched mean field theory and the partition 
function relations can be tested numerically. The easiest thing to check is the implication 
of eqn. (10) that the gravitational transition is independent of A. As in the case of pure 
gravity [4] , the expectation value of the number of nodes scales with the volume as 

< Af >~ N$ (44) 
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where 6 is one in the phase with finite coordination number and less than one in the 
phase where the coordination number diverges. For A = 2, we computed 5 for a = 1.6, 
a = 1.8, and a = 2.1 on lattices of volume 2000, 4000, 8000, and 16000. The results for 
5 are 0.55(1), 0.84(2), and 1.007(1) respectively. The error bars are statistical only. The 
confidence level of the (log) fits to a straight line were not high, suggesting that larger 
lattices are necessary to reach scaling, but the trend is consistent with a gravitational 
transition still in the region a = 1.91 — 2.16 of the A = transition. 

The gauge transition is more difficult. To search for the gauge transitions we computed 
the gauge part of the specific heat as a function of A for various fixed values of a on both 
sides of the gravitational transition and looked for a peak growing with volume. Figures 
5 and 6 show the results. There is a maximum in the specific heat in the region predicted 
by mean field theory, but there is no sharp peak, it does not grow with volume, and there 
is no evidence for critical slowing down. The same results were obtained for large and 
smaller values of a. This suggests that the gauge sector is in a strong coupling phase for 
all A. Other operators show consistent behavior. For instance, the expectation value of 
the plaquette is almost a straight line as a function of A in the region of A we considered 
and had no significant volume dependence on lattices up to a volume of 16000. One could 
attempt to verify that the entire phase is in strong coupling by examining Wilson loops 
for large A and testing for an area law, but this is much more difficult on a dynamical 
triangulation than on a fixed hypercubic lattice and we have not tried to do it. 

A simpler way to try to get information about what is going on in the gauge sector is 
to look at monopoles and their clusters. Consider a three-dimensional elementary volume 
(a tetrahedron here) with a surface of oriented plaquettes, Pi (i.e. P{ = 0). For each Pj 
define such that 

P t = 27m, + <t>i (45) 
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where — n < <j) < n. Then the monopole number in this volume is 



(46) 



Generally, a certain value of requires pi > (2nj — l)ir. Summing this inequality over the 
s sides of the elementary volume gives a limit on the monopole number: 



For a cube, this gives the familiar result that m E {—2,-1,0,1,2}. For a tetrahedron 
there are fewer possibilities: m G {—1, 0, 1}. A monopole density is defined as 



where N m is the number of tetrahedra with monopole number m and V is the volume (the 
number of tetrahedra). Clusters are defined as groups of neighboring monopoles (signs are 
ignored and monopoles are neighbors if the surfaces of their tetrahedra share a face). An 
interesting order parameter for clusters is the ratio 



where A max is the number of monopoles in the largest cluster and A^ot is the total number 
of monopoles. M is a good order parameter for a percolation transition, a transition from 
a phase with many clusters to a phase with one cluster. In the three-dimensional quantum 
gravity plus U(l) gauge fields system, the monopole density is 1/3 when A = and it 
slowly and smoothly decreases as A is increased. (On a cubic lattice the monopole density 
is 0.43 at A = and decreases for larger A.) The density has very little sensitivity to 
the value of a. The cluster order parameter, M, is small for A = and a in the region 
of the gravitational transition and it decreases smoothly as A is increased. Eventually it 
increases again, but this is only because the density becomes so small that there are of 
order unity monopoles (i.e. A max ~ -Ntot)- M also has no sensitivity to the value of a 



m < s/2 



(47) 



p= (N+ + N-)/V 



(48) 




(49) 
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and this is surprising. In general, one expects a low percolation threshold for lattices with 
a high coordination number (such as those generated in the negative curvature phase). 
Apparently, it is the coordination number of the tetrahedron that counts (since that is 
where the monopole lives) rather than that of the sites. A tetrahedron has a coordination 
number of four, the same as a site in a two-dimensional square lattice, so one might 
guess that the threshold for the tetrahedral lattice is roughly that of the square lattice, 
namely 0.59 [10]. This guess can be confirmed by artificially placing monopoles of a given 
density on the dynamical lattices and measuring the threshold directly. Consequently, the 
physical monopole density is always below the percolation threshold and they provide no 
useful signal. 

In conclusion, we have shown how to couple gauge fields to three-dimensional quantum 
gravity and made an effort to find an interesting continuum limit. No appropriate phase 
transition was found in the space of actions we considered. 

This work was supported, in part, by NSF grant PHY 92-00148. Some calculations 
were performed on the Florida State University Cray Y-MP and at the Pittsburgh Super- 
computer Center. We thank Olle Heinonen for discussions. 
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FIGURE CAPTIONS 



[1] The quenched mean field diagram for three-dimensional quantum gravity coupled 
to U(l) gauge fields. 

[2] The tetrahedra involved in a (1,4) move. 

[3] The tetrahedra involved in a (2,3) move. 

[4] The expectation value of the number of nodes, iVO, versus the ratio of the number 
of gauge update sweeps to the number of gravity update sweeps, a = 2.1, A = 1.0, and 
N3 = 4000. 

[5] The gauge component of the specific heat as a function of A at a = 1.8 for two 
lattice sizes: iV 3 = 2000 (crosses) and iV 3 = 4000 (boxes) . 

[6] The gauge component of the specific heat as a function of A at a = 2.2 for two 
lattice sizes: N3 = 2000 (crosses) and N3 = 4000 (boxes). 
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